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Computer  simulations  of  microwave  coupling  to  a  nitrogen-breakdown  plasma  have  been 
performed  at  25  Torr.  Nonhydrodynamic  ionization  fronts  are  observed  to  propagate  toward  the 
radiation  source  under  a  variety  of  circumstances.  Free-nitrogen-breakdown  simulations  in  a 
spherical  system  show  the  propagation  velocity  of  the  breakdown  wave  can  be  as  high  as  5  X 106 
cm/sec.  An  elementary  theory  is  used  for  estimating  the  speed  of  the  breakdown  wave  in  one 
dimension.  The  results  are  in  reasonable  agreement  with  breakdown  experiments. 
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PACS  numbers:  52.40.Db,  52.50.Dg,  52.80.Pi 


I.  INTRODUCTION 

The  pulsed  breakdown  in  air  and  other  gaseous  ele¬ 
ments  has  been  studied  extensively1  with  emphasis  on  the 
threshold  power  for  breakdown  and  its  dependence  on  the 
gas  pressure  and  radiation  wavelength.  The  hydrodynamic 
effects  were  first  considered  by  Lin  and  Theofilos.2  Their 
calculation  of  the  gas  heating  showed  that  both  high-field 
intensity  above  the  breakdown  intensity  and  increasing  mo¬ 
lecular  density  discouraged  energy  deposition  in  air.  Though 
the  preliminary  experimental  results  confirmed  the  general 
features  of  the  simple  theory,  the  observed  pressure  waves 
were  stronger  than  predicted.  Scharfman  et  al.3  have  ob¬ 
served  that  a  breakdown  plasma  prevents  the  transmission 
of  microwave  energy  to  points  beyond  the  plasma.  Experi¬ 
ments4  recently  performed  at  the  Naval  Research  Labora¬ 
tory  (NRL)  showed  that  the  microwaves  are  rapidly  decou¬ 
pled  from  a  target  or  plasma  surface.  This  rapid  decoupling 
of  the  microwaves  would  limit  the  amount  of  energy  that  can 
be  deposited  into  the  gas.  Computer  simulations  of  the  NRL 
experiments  have  been  performed  at  a  pressure  of  25  Torr  in 
nitrogen.  Throughout  this  work,  we  have  used  the  MINI 
code  developed  at  NRL  to  help  in  understanding  the  basic 
interaction  processes  in  the  experiments. 

II.  THE  MICROWAVE  NITROGEN  INTERACTION  CODE 
(MINI) 

The  microwave  nitrogen  interaction  code  (MINI)  used 
in  our  studies  is  a  one-dimensional  multispecies,  multitem¬ 
perature,  hydrodynamic,  nitrogen  chemistry  and  wave  op¬ 
tics  code.  The  species  followed  dynamically  were  N2,  N24 
(X ),  N24  ( B ),  N4f ,  N,  N(2D),  N2(/4  ’2),  N,(5  V),  N2(CV),  and 
the  electron  density  («,, ).  The  model  calculates  the  electron 
(7J,  vibrational  ( T„ ),  and  gas  (7^)  temperatures,  and  is  space 
and  time  dependent  with  microwave  absorption  and  reflec¬ 
tion  considered  in  the  wave-optics  mode.  The  detail  of  the 
chemistry  and  the  wave-optics  aspect  of  the  code  are  de¬ 
scribed  elsewhere*  and  will  not  be  repeated  here.  The  elec¬ 
trons  are  described  hydrodynamically  by  the  equations 

<?,(/»,)  +  <?y(fl,n,.)  =  Cr,  (1) 

ue  =  u„  —  prE  -  Dcdy( In  ne)  -  Dtdy( In  Te),  (2) 
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<3,10  +  a,  |  etut  +  p,  u,  +  Q,  j  =  R,-1I,  +  J,-E  +  Et 
f  t  =  Yn‘Tt  +  y  (3) 

where pe  =  neTe  is  the  electron  pressure,  je  —  —  eneut  is 
the  electron  current  density,  Rr  =  —  mtntvc(u,  —  um)  is 
the  electron  momentum  transfer  rate,  fie  =  e/mtvc  is  the 
electron  mobility,  and  De  =  Te/mtvc  is  the  electron  diffu¬ 
sion  coefficient.  The  electron  momentum  transfer  frequency 
is  vc ,  and  m,  is  the  mass  of  the  electron.  The  electron  energy 
flux  Q„  =  —  Krd/dy  Tt ,  where  the  electron  thermal  conduc- 
tivityis/r(  =fineTt  /mevc,  and/9  is  a  constant.  The  terms  on 
the  right-hand  side  of  the  energy  equation  represent  energy 
lost  due  to  heat  flow,  work  due  to  momentum  transfer,  joule 
heating  and  heating  due  to  elastic,  inelastic,  and  chemical 
processes.  The  term  Ce  represents  sources  and  sinks  terms 
for  the  electron  due  to  ionization,  recombination,  and  at¬ 
tachment.  The  Ee  term  in  the  electron-energy  equation  in¬ 
cludes  energy  losses  due  to  vibrational  and  electronic  excita¬ 
tion  of  molecular  nitrogen.  The  C,  and  Et  terms  are 
discussed  more  fully  in  Ref.  5. 

The  continuity  equation  for  each  of  the  heavy  particle 
species  is 

<?,  («,)  +  <?„(»,  «,)  =  C.-  (4) 

The  momentum  equation  for  the  heavy  particles  is 

d,(pu)  +  dy{pu2  +  PH  -(y/u  +i«b)<3,«I  =p+E-R„ 

(5) 

where  PH  =  In,  7^  is  the  heavy  particle  pressure,  p  =  2m, 
n,  is  the  mass  density,  p 4  =  2Z,e«,  is  the  positive  charge 
density,  and  p\pB)  is  the  kinematic  (bulk)  viscosity  coeffi¬ 
cients.  The  transport  coefficients  for  the  viscosity  and  ther¬ 
mal  conductivity  are  taken  from  Ref.  6.  The  heavy-particle 
energy  equation  is 

d,(eH)  +  dy\eHu  +  \P„  —  ( -j-A*  +/*s)<3, 1"  +  C?„  I  (6) 
=  nmT-,[e„(7;)-e„(7;)!  -jyE-Rcue  +EH 

=XC'"'7’*  +  Tpu2' 

i  *• 

where  eH  is  the  total  energy  density,  c,  is  the  translation- 
rotational  specific  heat,  £r(7’„)  =  ejexple/T',  )  —  1  j  “ 1  is  the 
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average  vibrational  energy  per  molecule,7  e~0.3  eV  is  the 
quantum  of  vibrational  energy  for  nitrogen,  r„  is  the  charac¬ 
teristic  time  for  vibrational  relaxation,8  and  E„  is  the  heat¬ 
ing  terms  due  to  chemical,  elastic,  and  inelastic  processes. 
Again  a  more  detailed  discussion  of  E„  can  be  found  in  Ref. 
'5.*Th£  heavy-particle  energy  flux  is  Q„  =  —  KHd/dy  Tg ,  and 
kh  is  theT&Epngl  conductivity  of  the  gas,  including  transla¬ 
tional  and  roqponal  contribution.9  The  mass-average-ve¬ 
locity  of  ththeayy  particles  is 

«=  —  ~um. 

P  > 

The  vibrational  energy  equation  may  be  written  ap¬ 
proximately  as 

d,[nme„(Tv)\  +  dy\nmev(Tv)um  +Qn] 

=  -  ejr,  ))+£„, 

where  Qv=  —  Kudy  T„  is  the  vibrational  energy  flux  and  kv 
is  the  vibrational  conductivity.  The  vibrational  temperature 
can  be  calculated  assuming  that  the  nitrogen  molecule  is  a 
harmonic  oscillator  and  that  the  vibrational  levels  have  a 
Boltzman  distribution.  Hence,  the  vibrational  energy  source 
term  due  to  electron  impact  is10 

E„  =  enmnt  { 1  -  exp(  -  e/T.)] 

2S_1«*.|1-  exPMr„  -  Te)/T„  T,  j },  (7) 

where  Xv  is  the  excitation  rate  coefficient10  for  the  irth  vibra¬ 
tional  level  obtained  from  the  experimentally  measured 
cross  sections. 11 

The  dc  portion  of  the  electric  field  and  the  diffusion 
velocities  are  determined  by  demanding  the  total  current 

Jr  =  J,  +  X  =  0*  (8) 

i 

where  J,  =  Z,entu,  is  the  current  density  of  the  ith  ion  spe¬ 
cies.  In  steady  state,  each  of  the  ion  species  satisfies  an  equa¬ 
tion: 

u,  =  um  +  /r,E  -  D,dy{ In  «,)  -  D,dy[ In  TJ,  (9) 

where  the  ion  diffusion  coefficient  and  mobility  are  related  to 
the  collision  frequency  v,  by  D,  =  Tg/M,vt  and  y,  =  Z,e/ 
nt/V,.  The  mass  of  the  /th  ion  species  is  m,  .  Summing  Eq.  (9) 
for  all  ion  species  and  using  Eq.  (8)  gives  an  equation  for  the 
ambipolar  electric  field.  Once  the  ambipolar  field  is  found, 
Etfc  (9)  and  (2)  are  used  to  solve  for  the  flow  velocities  of  the 
ions  id  the  electron. 

'the  Eulerian  difference  equations  of  Eqs.  (1)— (7)  are 
solved  using  the  algorithm  of  Rubin  and  Burstein.12  The 
scheme  is  simple,  second  order  in  accuracy,  and  is  relatively 
stable  while  centering  the  dissipation  terms.  The  hydrodyna¬ 
mics  are  determined  from  the  updated  values  of  the  chemis¬ 
try.  Hence,  the  hydrodynamics  are  considered  as  a  small 
correction  to  the  chemistry. 


II.  CASCADE  IONIZATION  OF  NITROGEN  BY  A 
MICROWAVE  PULSE 

Microwave  gas  breakdown  begins  when  a  small  number 
of  priming  electrons  acquire  sufficient  energy  from  the  elec¬ 


tric  field  to  ionize  the  gas.  The  priming  electrons  gain  energy 
when  the  oscillatory  energy  of  the  electrons  in  the  electric 
field  is  randomized  by  collisions  into  thermal  energy.  Once 
an  electron  has  an  energy  in  excess  of  the  ionization  energy  of 
the  gas,  the  electron  can  easily  ionize  the  gas,  resulting  in  the 
generation  of  two  low-energy  electrons.  The  electron  ava¬ 
lanche  can  be  described  by 

*«(0  =  n«(0)exp(r/rJI),  (10) 

where  ne{t)  is  the  electron  density  at  time  t,  and  rB  is  the 
breakdown  time.  Physically,  rB  is  the  time  required  for  an 
electron  to  gain  the  energy  needed  to  ionize  the  gas.13  Figure 
1  shows  the  breakdown  time  tb  (sec)  and  the  electron  tem¬ 
perature  r,(eV)  for  conditions  typical  of  microwave  break¬ 
down  of  Nz.  A  high  radiation  field  would  promote  the  elec¬ 
tron  avalanche  by  increasing  the  energy  imparted  to  the 
electron  per  collision  with  the  molecule.  A  high  molecular 
density,  however,  requires  higher  intensities  for  the  avalance 
breakdown.  At  moderate  pressures  above  25  Torr,  the  domi¬ 
nant  energy  loss  mechanism  for  the  electrons  is  the  excita¬ 
tion  of  the  electronic  and  vibrational  modes  of  the  gas.  The 
root-mean-square  electric  field  to  density  ratio  E/nm 
(V  cm2)  does  not  uniquely  define  the  absorbed  energy  for  an 
ac  electric  field.  The  striped  region  in  Fig.  1  is  bounded  by 
the  constant  intensity  line  (I)  and  the  constant  pressure  line 
(P).  The  rapid  variation  in  tb  with  increasing  intensity  allows 
a  breakdown  threshold,  IB,  to  be  assigned  at  the  onset  of 
breakdown.  The  ability  to  couple  microwave  energy  to  a 
breakdown  gas  depends  sensitively  on  the  breakdown  time, 
tb,  and  the  electron  density. 


FIG.  1 .  The  breakdown  time  and  the  electron  temperature  vs  the  rms  elec¬ 
tric  field  to  density  ratio  E  /n„ . 
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The  efficiency  of  converting  microwave  energy  to  the 
translational  energy  of  the  gas  can  be  enhanced  by  coupling 
the  microwaves  to  a  preformed  plasma.  Plasma  maintenance 
of  the  preformed  plasma  by  microwaves  allows  the  micro- 
wave  deposition  to  be  localized  while  maintaining  a  high 
electron  density.  In  short-pulse  breakdown  experiments, 
with  pulse  lengths  tp  <  1  //sec,  a  high  electron  density  is  a 
necessary  condition  for  rapid  heating  of  the  gas.  Plasma 
maintenance  requires  the  breakdown  time  rB  >  rp  when  ne 
(0)>  1 .  If  the  local  electric  field  is  not  carefully  matched  to  the 
instantaneous  plasma  conditions,  the  gas  heating  will  be  self¬ 
limited  by  the  formation  of  a  nonhydrodynamic  ionization 
front.  Gas  breakdown  off  a  reflecting  surface  without  a  pre¬ 
formed  plasma  will  not  result  in  a  local  deposition  of  the 
microwave  energy.  Since  the  microwave  deposition  is  deter¬ 
mined  self-consistently  with  the  electron  density  profile,  the 
resulting  plasma  density  cannot  rapidly  heat  the  gas  in  the 
short-pulse  experiments. 

A.  Gas  breakdown  near  a  reflecting  surface 

A  plane  wave  incident  normal  to  a  reflecting  surface 
requires  less  power  for  gas  breakdown.  The  constructive  in¬ 
terference  of  the  incident  and  reflected  microwave  results  in 
the  characteristic  standing  wave  pattern  for  the  electric  field 
and  a  factor  of  4  increase  in  the  field  intensity.  If  the  intensity 
is  above  the  breakdown  threshold,  breakdown  of  the  gas  oc¬ 
curs  intantaneously  at  the  peak  nodes  of  the  electric  field.  As 
the  electron  density  rises,  these  “pancakes”  of  electrons  be¬ 
gin  to  attenuate  the  radiation.  Pancakes  closest  to  the  reflect¬ 
ing  surface  rise  more  rapidly  early  in  the  pulse  (Fig.  2).  As  the 


electron  density  continues  to  rise,  attenuation  of  the  micro- 
wave  causes  the  pancakes  closest  to  the  surface  to  fall  in 
density  below  that  of  pancakes  far  from  the  surface.  Shorter 
systems  with  fewer  pancakes  decouple  more  rapidly  than 
longer  systems.  However,  the  total  absorbed  energy  remains 
essentially  constant,  independent  of  the  system  size  for  a 
fixed  pressure  and  pulse  length.14 

The  finite-length  simulations  represent  the  planar  re¬ 
gion  near  the  surface  of  a  focused  microwave  system  where 
the  distance  from  the  surface  is  less  than  the  radius  of  the 
spotsize.  Figure  2  shows  two  similar  simulations  with 
I  =  6.25  kW/cm2  for  a  one-and-a-half  and  three  wavelength 
system  without  hydrodynamic  effects.  The  plasma  absorbs 
approximately  20%  of  the  total  microwave  energy  with  an 
asymptotic  absorption  efficiency  of  80%.  The  absorption  is 
regulated  by  the  electron  density.  Shorter  systems  have  a 
much  higher  electron  density  than  larger  systems  and  can 
more  readily  heat  the  gas  in  the  short-pulse  experiments. 
The  gas  is  not  expected  to  be  heated  more  than  the  observed 
4%  above  the  ambient  temperature  in  a  weakly  focused  sys¬ 
tem.  The  absorbed  energy  is  distributed  over  a  larger  number 
of  pancakes  with  a  very  low  electron  density.  As  will  be  dis¬ 
cussed  later,  the  low  electron  density  of  Fig.  2  results  in  a 
long  heating  time  for  the  gas  of  rg  ~6.5  fi  sec,  where  rg = Tg 
(dTg  /dt)~'.  The  time  required  to  heat  the  gas  to  an  apprecia¬ 
ble  temperature  in  a  weakly  focused  system  would  require 
pulse  lengths  ofrp>6.5//  sec.  However,  even  for  long  pulses, 
a  significant  heating  of  the  gas  is  never  realized.  Later  in  the 
pulse,  the  pancakes  furthest  away  from  the  reflecting  surface 
continue  to  rise  in  density.  Ultimately,  a  single  pancake  will 
have  sufficient  density  to  reflect  the  radiation.  This  pancake 


1280 


J.  App).  Phya.,  VcH.  54,  No.  3.  March  1983 


Yee,  Ali,  and  Bollen 


1280 


acts  as  a  secondary  reflecting  surface  and  decouples  the  radi¬ 
ation  from  the  primary  metallic  surface.  Gas-breakdown  ex¬ 
periments4  using  a  reflecting  surface  show  this  behavior. 
Similarly,  free-nitrogen-breakdown  simulations  without  a 
surface,  using  a  focused  system,  also  show  that  the  radiation 
is  decoupled  from  the  initial  breakdown  site. 

B.  Maintenance  of  a  preformed  plasma  by  microwaves 

Microwave  coupling  to  a  preionized  plasma  near  a  me¬ 
tallic  surface  show  some  general  characteristics.  Though  the 
plasma  absorbs  the  microwave  energy,  the  absorption  can¬ 
not  be  localized.  A  traveling  ionization  front  is  observed  to 
proceed  toward  the  radiation  source.  Figure  3  shows  the 
time  history  of  the  ionization  front  in  the  strong  mismatch 
case  7>/|,  and  rp  >  rc.  The  critical  time  rc  is  the  time  re¬ 
quired  for  the  electron  avalanche  to  reach  the  critical  density 
nc^(o2me/Aire2c^  10I3/A  2(cm),  where /l  is  the  wavelength  of 
the  incident  radiation.  In  this  simulation  I  =  12.5  kW/ 
cmV,  =  3  kW/cm2,  rp  =  0.5  //sec,  rc  =  0. 16  //sec,  and 
all  hydrodynamic  effects  are  suppressed.  The  initial  electron 
density  (preionized  plasma)  was  initialized  with  a  symmetric 
profile  with  a  peak  density  of  0. 1  nc  and  full  width  half¬ 
maximum  ofO.lA.In  Fig.  3,  the  peak  density  is  determined 
by  Eq.  ( 10).  However,  the  breakdown  time  tb  is  a  complicat¬ 
ed  function  of  the  incident  and  reflected  microwaves.  It  is 
clear  that  the  small  reflectivity  (I?  <  10%)  is  important  in 
determining  the  propagation  speed  since,  as  Fig.  3  shows,  the 
separation  distance  between  peaks  is  A  /4.  The  ionization 
front  proceeds  toward  the  source  at  a  characteristic  speed 
D  (t )  =  A  /4rc(r ).  The  maximum  density  is  decreasing  with 
time  while  the  total  absorption  continues  to  increase  with 
time.  The  asymptotic  efficiency  is  approximately  95%.  The 


movement  of  the  ionization  front  limits  the  time  the  gas- 
plasma  system  samples  the  radiation  field.  The  sampling 
time  rs  is  determined  by  a(cm)  =  SlsD  {t )  dt~D  (0)rs,  where 
D  (0)  =  A  /4tc(0)~1.6x  106  cm/sec  is  the  early-breakdown 
speed.  Estimating  the  absorption  length  as  a(cm)~  A  /2ir 
gives  a  sampling  time  of  rs~rc.  In  the  pressure  regime 
P>  25  Torr,  the  gas  is  primarly  heated  through  the  quench¬ 
ing  of  the  electronic  states  by  molecular  and  atomic  nitro¬ 
gen.  The  fractional  change  in  the  gas  temperature  is  eg 
=^TS/TS) =  Ts2afiQa?nan0€a/cunm  Ts,  where  Qaf,  is 
the  quenching  rate  of  the  N2  triplet  state  na  by  np.  The  spe¬ 
cific  heat  of  nitrogen  is  cu ,  and  ea  is  the  energy  level  of  the  na 
state.  As  an  example,  the  N2(CV)  state  is  quenched  by  N2 
with  a  rate  coefficient15  of  1.2x10“"  cmVsec.  The  frac¬ 
tional  change  in  the  gas  temperature  eg  is  approximately 
21%  from  the  ambient  temperature. 

The  region  of  the  maximum  gas  heating  is  at  the  initial 
breakdown  site  ( ~  A  /4  from  the  surface).  The  electric  field  is 
less  than  1  %  of  the  incident  time-average  electric  field  at  the 
maximum  gas  temperature  [Fig.  4(a)].  The  plasma  absorbs 
approximately  70%  of  the  incident  microwave  energy  and 
has  an  electron  temperature  of  less  than  2.5  eV  [Fig.  4(b)]. 
Maximum  gas  heating  requires  the  sampling  time  and  the 
electron  temperature  to  be  as  large  as  possible  without  de¬ 
coupling  the  system.  A  high  electron  temperature  populates 
the  electronic  states  which  are  rapidly  quenched  to  heat  the 
gas.  Approximately  4%  and  40%  of  the  incident  energy  goes 
into  heating  the  gas  and  exciting  the  vibrational  states.  A 
similar  simulation  with  the  peak  power  reduced  to  1  =  6.25 
kW/cm2  shows  no  movement  of  the  initial  ionization  region 
while  increasing  the  fractional  change  in  the  gas  temperature 
to  27%.  Further  irradiance  would  increase  the  gas  tempera¬ 
ture. 
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FIG.  3.  Time  history  of  the  ionization 
front  in  a  planar  system  showing  the  wave 
motion  toward  the  radiation  source.  The 
microwaves  are  incident  from  the  left  onto 
a  preionized  plasma  near  a  reflecting  sur¬ 
face  on  the  right  boundary.  The  normali¬ 
zation  parameters  are  the  same  as  in  Fig. 
2. 
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FIG.  4.  (a)  The  time-averaged  electric  field  square  and  (b)  the  electron  tem¬ 
perature  vs  position.  The  normalization  parameters  are:  TA  —  0.25  eV, 
*  =  7.33  cm  - and  the  arbitrary  units  of  <£2>  =  1  correspond  to  an  aver¬ 
age  intensity  of  9.4  kW/cm2. 


C.  Gas  breakdown  in  a  focused  microwave  system 

Simulations  of  a  focused  microwave  system  were  also 
performed.  If  the  field  intensity  in  a  focused  sys*em  greatly 
exceeds  the  breakdown  threshold  at  the  focal  spot,  a  break¬ 
down  wave  propagates  toward  the  source.16  The  motion  of 
the  microwave  plasma  absorbing  layer  is  not  a  simple-break¬ 
down  wave,  since  the  reflected  microwave  is  important  in 
determining  the  propagation  speed.  Figure  5  shows  the  time 
history  of  the  electron  density  profile  during  the  course  of  a 


rp  =  0.5-^sec  pulse  in  a  focused  system.  The  spherical  wave 
equation  is  solved  for  the  electric  field  to  model  the  conver¬ 
gence  of  the  field.  The  peak  incident  field  is  7=  10  kW/cm2 
with  a  system  aspect  ratio  of  2.5.  The  average  power  at  focus 
is  31.25  kW/cm2,  and  all  hydrodynamic  effects  are  sup¬ 
pressed.  The  very  rapid  early-time  speed  of  the  absorbing 
wave,  Z>~5X  106  cm/sec,  gives  a  sampling  time  of  rs 
~1.7x  10-8  sec.  However,  the  gas  continues  to  be  heated 
after  the  passage  of  the  field  by  the  de-excitation  of  the  elec¬ 
tronic  states.  The  miximum  gas  and  vibrational  temperature 
are  e,~10%  and  e„~260%,  respectively.  The  quarter- 
wavelength  fine  structure  in  the  density  profile  shows  the 
effected  light.  The  time-average  reflection  and  absorption 
coefficient  are  19%  and  68%,  respectively.  A  snapshot  of 
the  electric  field  and  electron  temperature  (Fig.  6)  shows  the 
lack  of  field  penetration  (ar~0. 1 A )  and  the  cold  plasma  tem¬ 
perature,  Tr  <  3.3  eV.  The  gas  and  vibration  temperatures 
again  are  reduced  when  hydrodynamic  effects  are  included 
in  the  simulations.  The  gas  is  expected  to  be  heated,  at  most, 
to  10%  of  the  initial  value  during  the  course  of  a  0.5-/xsec 
pulse. 

The  propagation  velocity  of  the  breakdown  wave  in  a 
focused  microwave  system  can  be  calculated.  Simulations 
show  the  absorption  is  localized  at  the  head  of  the  break¬ 
down  wave.  The  reflection  coefficient  is  relatively  constant 
after  the  initial  breakdown  (t  >  rc).  The  breakdown  time  rB 
is  inversely  proportional  to  the  light  flux  l/rB  oc 
(1  +  2  Re£  +  £  V,  where  the  reflection  coefficient  is£  2.  The 
intensity  in  a  conically  focused  system  can  be  written  as 


/(*,»)« 


d  +  2iy+ttxm 


where  ^ (f )  is  a  form  function  characterizing  the  pulse.  The 
radius  of  the  radiation  channel  at  position  X  on  the  channel 
axis  is  r  =  r0  4-  X  tan  0.  The  minimum  radius  of  the  waist  at 
the  focal  point  is  r0,  and  6  is  the  half-angle  of  the  radiation 
column.  Following  Raizer,16  the  propagation  velocity  of  the 


FIG.  5.  Time  history  of  the  ion¬ 
ization  front  in  a  spherical  sys¬ 
tem.  '’’he  microwaves  are  inci¬ 
dent  from  the  left  boundary 
with  an  intensity  of  5  kW/  cm2. 
The  focused  average  power  is 
3 1 .25  kW/cm2  at  the  right-side 
boundary. 
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FIG.  6.  (a)  The  time-averaged  electric  field  square  and  (b)  the  electron  tem¬ 
perature  vs  position  for  the  simulation  of  Fig.  5.  The  normalization  param¬ 
eters  are  the  same  as  in  Fig.  4. 


breakdown  wave  is  approximately  D  =  r0/rc  tan  6.  Under 
the  assumptions  made,  the  velocity  D  would  be  an  upper 
bound  on  the  propagation  velocity  of  the  breakdown  wave. 

IV.  SIMULATION  OF  THE  NRL  EXPERIMENTS 

We  compare  the  simulation  results  with  the  NRL  ex¬ 
periment.  The  NRL  35-GHz  1 12-kW  focused  gyro'  .  _ys- 
tem  gives  an  average  power  at  the  focal  spot  of  33  kW/cm2. 
The  polarization  of  the  microwave  is  such  that  the  electric 
field  is  perpendicular  to  the  plasma  gradient  (S  polarization). 
The  average  aspect  ratio  of  the  experiment  gives  an  equiva¬ 
lent  simulation  aspect  ratio  of  2.5.  The  free-nitrogen  toreak- 
down  experiments  were  performed  at  25  Torr  with  a  100  pps 
repetition  rate  at  rp  =  1.5/xsec.  The  minimum  radius  of  the 
waist  at  the  focal  spot  is  r0  —  0.75  cm,  and  the  half-angle  of 
the  radiation  column  satisfies  tan  9  =  0.35.  Experimentally, 
framing  camera  data  (Fig.  7)  gives  a  propagation  velocity  of 
Z>~4.6X  10*  cm/sec.  The  width  of  the  absorbing  wave  is 
approximately  1.2  cm.  The  measured  absolute  intensity  of 
the  second  positive  band,  N?  2 P  (0-0)  at  337 1  A,  is  1 .4  X  1 02 1 
eV/cm2  sec.  The  calculated  spatial  width  of  the  337 1  -  A  light 
emission  is  0.8  cm,  with  a  peak  intensity  of  1 .8  X 102 1  eV/cm2 
sec.  As  mentioned  previously,  the  calculated  speed  of  the 
breakdown  wave  is  5  X 106  cm/sec.  Experimentally,  the  in¬ 
ferred  breakdown  time  is  ra~2.5x  10~"  sec  or  rc 
~7.4x  10~7  sec.  This  is  to  be  compared  to  the  calculated 
breakdown  time  early  in  the  pulse,  ra~4.5x  \0_<)  sec,  and 
late  in  the  pulse,  rB^ 6X  10"*  sec.  The  simulation  results 
are  felt  to  be  in  reasonable  agreement  with  the  experiment. 


FIG.  7.  Framing  camera  data  showing 
the  motion  of  the  ionization  front  in  a 
focused  free-nitrogen-breakdown  ex¬ 
periment.  The  top  and  bottom  photo¬ 
graphs  were  taken  0.2  and  0.1  fisec 
after  the  first  observable 
breakdown. 


I - 1 

10  cm 

V.  CONCLUSION 

In  conclusion,  computer  simulation  of  microwave  cou¬ 
pling  to  a  weakly  ionized  gas  plasma  shows  excellent  absorp¬ 
tion  of  the  microwave  energy.  The  electronic  states  of  nitro¬ 
gen  excited  by  electron  impact  are  rapidly  quenched  by 
molecular  and  atomic  nitrogen  to  heat  the  gas  at  pressure 
above  25  Torr.  However,  high  power  irradiance  can  easily 
produce  a  non-hydrodynamic  ionization  front  which  can  se¬ 
verely  limit  the  heating  of  the  gas.  Microwave  gas  break¬ 
down  experiments  with  pulse  lengths  less  than  one  micro¬ 
second  are  expected  to  show  only  a  slight  heating  of  the  gas. 
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